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We studied the edge states and transverse electron focusing in the presence of spin-orbit interaction 
in a two dimensional electron gas. Assuming strong spin-orbit coupling we derived semiclassical 
quantization conditions to describe the dispersion of the edge states. Using the disprsion relation 
we then make predictions about certain properties of the focusing spectrum. Comparison of our 
analytical results with quantum mechanical transport calculations reveals that certain features of 
the focusing spectrum can be quite well understood in terms of the interference of the edge states 
while the explanation of other features seems to require a different approach. 
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I. INTRODUCTION 

Semiclassical approximations are often very practi- 
cal for understanding certain physical phenomena. Be- 
yond their practicality, they also provide a rather general 
framework to treat quantum systems of interest. A good 
example to illustrate the merits of semiclassical approach 
is the transverse electron focusing (TEF). The geometry 
of electron focusing is shown in FigHJ The current is in- 
jected into the sample at a quantum point contact called 
injector (I) in perpendicular magnetic field B. If the mag- 
netic field is an integer multiple of a focusing field B focus , 
electrons injected within a small angle around the per- 
pendicular direction to the edge of the sample can be 
focused onto a collector quantum point contact (denoted 
by C in Fig. [T]) which acts as a voltage probe. Therefore, 
if the collector voltage is plotted as a function of mag- 
netic field one can observe equidistant peaks at magnetic 
fields B = p * B focus {p = 1,2,3,...) corresponding to 
cases where the cyclotron diameter 2Rc{B) is an integer 
multiple of the distance L between the injector and the 
collector. 

TEF is a versatile experimental technique (for a re- 
view of the various problems where it has been used see 
Ref. [ll). In the case of quantum wells containing two di- 
mensional electron gas (2DEG), the accessibility of the 
quantum ballistic transport regime opened up the way 
to the experimental demonstration of coherent electron 
focusing'^1^ as well. Recently, several experiments inves- 
tigated the effect of spin-split bands in semiconductors 
on magnetic focusing^"—. Of special interest are for us 
the experiments of Refs. in which evidence of spin- 
orbit interaction (SOI) dependent focusing have been 
found. These experiments sparked considerable theoret- 
ical interest^i^ d^d^ as well. Refs. Hiinilil have in com- 
mon that they consider the properties of bulk Landau 
levels in the presence of SOI to explain the experimen- 
tal results on focusing. While using the bulk Landau 
levels as a starting point is certainly justified when dis- 
cussing magneto-oscillation&i^^'*, for the geometry shown 
in Fig. [T] one expects that the edge states should play a 
central role in the transport phenomena. Indeed, this 




FIG. 1: (color online) Schematic geometry of the transverse 
electron focusing setup. The 2DEG is contacted by an injector 
(I) and a collector (C) probe and perpendicular magnetic field 
is applied. Classical quasiparticle trajectories leaving from the 
injector at normal direction, depending on the strenght of the 
magnetic field, can be focused onto the collector. The inset 
shows the details of the tight-binding model used to describe 
the the probes (c.f. the inset of Fig. 1 in Ref. 0). 



was the approach adopted in Ref. to discuss coher- 
ent electron focusing. The rich physics brought about 
by the interplay of SOI and the confinement due to ex- 
ternal magnetic field and electrostatic potential has also 
attracted significant theoretical attentioni^"— but impli- 
cations on electron focusing have not been considered. 

Here we aim to investigate whether the electron focus- 
ing spectrum in 2DEG with strong SOI can be explained 
in terms of edge states formed as a combined effect of 
SOI, magnetic field and (an assumed) hard wall confine- 
ment potential. To this end we first derive semiclassical 
quantization conditions which describe the dispersion re- 
lation of the edge states in the limit of strong SOI and 
weak magnetic fields. These results shed new light on and 
help to better understand the exact quantum solution of 
this problem, published very recently in Ref. H^. We then 
study how the properties of the edge states are mani- 
fested in the transport phenomena of the focusing setup 
shown in Fig.[T] We expect that our results should be rel- 
evant in the case of e.g. InSb quantum wells, where the- 
oretical predictions^i and recent experiments^>22ri^ indi- 
cate that it is possible to fabricate samples with strong 
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(compared to GaAs/AlGaAs heterostructures) spin orbit 
interaction and ballistic quasiparticle propagation over 
distances of the order of 1/xm at low temperatures. 

The rest of the paper is organized in the following way: 
in the next section we briefly introduce the semiclassical 
framework that we will be using. In Section IIIII we de- 
rive a pair of semiclassical quantization conditions for the 
edge states and compare the obtained band structure to 
the results of exact numerical calculations. We also dis- 
cuss how our semiclassical results are related to other ap- 
proximation methods found in the literature. The semi- 
classical quantization conditions then allows us in Sec- 
tion |IV] to make prediction about the focusing spectrum. 
We end our paper by a comparison of these predictions 
to numerical transport calculations and a short summary 
in Section IVl 



II. SEMICLASSICAL THEORY WITH SPIN 
DEGREES OF FREEDOM 

Generally, the non-relativistic single-particle Hamilto- 
nian of spin 1/2 particles can be written as^ 



H — Ho + Hi , Ho 



2m* 



V{v) 



(1) 



where r and p denotes the position and momentum op- 
erators, respectively, m* is the effective mass of the par- 
ticles and we assume that the spin dependent part have 
the following form: 



Hi = /iK C(r, p) • (T. 



(2) 



Here cr = [ax, o'yjO'z) is a vector of Pauli matrices and the 
h comes from the spin operator s = ^ha. The constant 
K gives the spin orbit coupling. 

The application of semiclassical methods developed for 
systems that can be described by scalar Hamiltonians^^ 
is not straightforward if one has to consider the spin de- 
gree of freedom as well. Namely, as a first step one would 
have to define a classical Hamiltonian which is not a triv- 
ial task since there is no classical analogue of the spin. 
Various semiclassical schemes have been proposed, see 
Refs. [ToH28l - [30 l and we refer to these original papers for 
most of the details. A short overview of the different 
approaches can also be found in Refs. ilia. 

Here we will follow the approach first used by Yabana 
and Horiuchi21 and later generalized and further devel- 
oped in Refs. [28l - l32l . It is often called the strong coupling 
limif^ because it corresponds to a double limit /i — )■ and 
K — >■ oo while k — Hk is kept constant^. For spin 1 /2 par- 
ticles this approximation scheme leads to two classical 
Hamiltonians 



=na±k\C{r,p) 



(3) 



where the vector C(r, p) is the phase-space symbol of the 
operator C(r, p)2&i22, and represents an effective magnetic 



field which depends on the classical variables r, p. This 
approximation introduces semiclassical phase corrections 
to the orbital motion^li^^. It is restricted however to the 
case of |C(r, p)| > because at phase-space points where 
|C(r,p)| vanishes and therefore V.^ becomes degenerate, 
mode conversion between trajectories described by Ti^ 
occurs, posing a serious difficulty to the theory. 

For a two dimensional electron gas (2DEG) in perpen- 
dicular magnetic field and assuming that Hi describes 
Rashba-type spin-orbit coupling, which will be our main 
interest in the rest of the paper, the "strong coupling" 
approach results in the following semiclassical Landau- 
level spectrum^!^: 



1,2, 



(4) 



where ujc — ^ is the classical cyclotron frequency, Ib = 
y/h/eB is the magnetic length, and by using the notation 
an for the coupling constant R in this particular case, kso 
is given by kg^ — m*aji/h. Comparing this to the exact 
result^^ 



Eo 

En 



n ± \ 2n{ksoh 



(5) 



we see that the semiclassical Landau levels are good ap- 
proximations of the exact ones if 2n{ksolB)'^ ^ 1/4. This 
requires large quantum numbers n (i.e. large Fermi en- 
ergy) and/or strong spin-orbit coupling and not too 
strong magnetic field (i.e. Ib ^ I/a/S is not too small). 
We expect therefore that the strong coupling method 
should be adequate if these conditions are met. An es- 
timate of the appropriate magnetic field range assuming 
InSb quantum well material parameters will be given af- 
ter Eq. (121]). 



III. SEMICLASSICAL THEORY OF EDGE 
STATES 

We assume that the 2DEG, formed e.g. in the quan- 
tum well of an InSb heterostructure, is in a perpendicular 
homogeneous magnetic field. The motion of electrons is 
confined by a hard- wall potential, V{y) = oo for y < 
(see Figll] for the geometry). The quantum mechani- 
cal description of the system can be obtained using the 
Hamiltonian (|Tj) where Ho — corresponds to the ki- 
netic energy of particles and the operator tt is defined as 

= (t^x, T^y) — p + eA, where p = —ih^ is the canonical 
momentum operator and A is the vector potential. Fur- 
thermore, Hi = an,{TTx(Ty — ■KyUx) describes the Rashba 
spin-orbit (RSO) coupling in the system, ax, <Jy are Pauli 
matrices acting in the spin space. We assume that an is 
constant in space and neglect its possible random varia- 
ton due to nanosize domainsSS. 

To preserve the translational invariance of the system, 
we choose the Landau gauge A = [By, 0, 0)"^. Using the 



3 



ansatz ^'(r) = e*'^^<i>(t/) we can simplify the problem to 
an effectively one dimensional (ID) one, which we will 
solve in semiclassical approximation. The discussion goes 
along the lines of Refs. 31 [for a recent application 

see also Refs. [37ll38| . We seek the solutions of the ID 
Schrodinger equation H^{y) = E^{y) in the following 
form—: 



E 

q>0 



(6) 



where aq{y) are spinors and S{y) is the classical ac- 
tion. Performing the unitary transformation $ — >■ 
e-HSto)$(y), H -> e-i^'^y^HeT^S<^y^ the Schrodinger 
equation can be rewritten as 



2m* 



E 

lily) 



(My) 



2m* 



(7) 



= 0. 



Here 11^ = 11° = + eA^ir), being = hk, and 
ily=Py+Wy, where Wy = ^|^. The WKB strategy^e 
is to satisfy Eq. ([7]) separately order by order in h. 
At 0{h^) order we obtain 



-^^-"^ ^n^oiSV''^^ Uo(.)-0. (8) 



E 



2m* 

Nontrivial zeroth order eigenvectors ao(2/) exist if 

2m* 



(9) 



where Q = ^(ng)2 + (n0)2. This means that Q is a 

constant of motion for a given energy E and the two 
branches of Eq. ^ define 



Q+ = pso + \Jvlo + 2m* 



(10a) 



(10b) 



where we used the notation p^^ — m*aji. We find there- 
fore that the classical equations of motion that can be 
derived from Eqs. ^ and (jlOp represent two harmonic 
oscillators. The corresponding zeroth order eigenvectors 
of Eq. dH) are 



V± = 



V2 ( 



62 



f(e(l;)T^/2) 
(0(y)T7r/2) 



(11) 



can be obtained from the 0(h}-) order of Eq. ([T]). Using 
Eq. ([TTT) we find that the 0(h}) order equation can be 
cast into the following form: 



no 

X 

no 

y 



cos 9{y) 
sm9{y) 



0. (12) 



This equation does not depend on 7± , which means that 
with the choice of the eigenvectors shown in Eq. pT|) , the 
phase 7±(j/) is already determined up to an unimportant 
constant factor. Moreover, by rewriting Eq. as 



(13) 



where V,^ — /2m? ± aji Q, it is easy to see that it 
expresses probability current conservation and it can be 
solved for A± using standard methods^^. From these 
results one finds that $(?/) in semiclassical approximation 
is given by 



1 



2in 



0,±| 



e5(e*fe)TT/2)) 



(14) 



where \^y'^ 



(n0)2. nowever, as the momen- 
tum tty'^ is a multi-valued function, we need to intro- 
duce the index j = +1,-1 do distinguish the different 
branches. The corresponding classical actions S^{y,y^) 
read 

Sfiy, y^) = 3 r ^Ql-miyTdy' (15) 

•'Vo 

where as usually, we have chosen the classical turning 
points as the phase reference points for the action. Sim- 
ilarly, the phase 0^{y) = &fiy) is multivalued as well. 

We now have to take into account the confinement po- 
tential V{y). The transverse wavefunctions ^{y) shown 
in Eq. (jl4p would not satisfy the boundary condition at 
y = 0. In order that the transverse wavef unction does 
satisfy the boundary condition we make a linear com- 
bination of the functions ^f{y) defined above. We try 
the following ansatz for the transverse semiclassical wave- 
function: 



My) 



e 4 cos 



|nr(y)| V"^*^ 





h 


2 




h ^ 


2 




(y-Va) 

h 


e-{y) 

2 


ho 


(y-yo) 1 






h 


2 



(16) 



where 9{y) is the phase of ^^{y) — i\i^{y). nowever, 
the eigenspinor ao(y)^ can more generally be sought as 
= ^±(y)e*'''±(^^V± where A±(y) is a real amplitude 
and 7±(j/) is a phase. Equations for A±(ti) and ^±{y) pear due to the magnetic field. The turning points are 



where C+, C_ are constants and the 7r/4 factor in the ar- 
gument of the cosine functions takes into account the 
effect of the classical turning points at y^ which ap- 
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given by the physically acceptable zeros of the equation 
- [n°(?/')]^ = 0. Note that |)(y) also depends on k 
through and 9^ but in order to keep the notations 
uncluttered we did not write this out explicitly. The dis- 
persion relation for the edge states can be obtained by 
demanding that the wave function vanishes at the bound- 
ary, i.e. = 0) = 0. This is a homogeneous system 
of equations and nontrivial solutions can only be found 
if the respective determinant is zero. This results in the 
following implicit dispersion relation: 



cos 



_ 6)+(0) 
2 

e'(o) , 



cos 



f ) COS ( 



h 



f q2i + 

e+(o) , 21 

2 ^4 



(17) 

Hence, if we denote by $'^(?/) and $ (y) the first and 
second spinors appearing in Eq. ([TO)) , respectively, the 
transverse wave function can be written, apart from a 
normalization factor, as ^ ^^(y) ^Tso ^ {y) where 
rso is given by 



cos 



e+(o) 



n?;+(0)|eos(^^%^-q^^ + f) 



(18) 

For fc w the reflection amplitude is ~ — e"/^. 

Furthermore, with the help of trigonometric identities 
and assuming that 



cos 



6)+(0) 0-(O) 



1 



(19) 

we find from Eq. (1171) the following quantization condi- 
tion: 



cos 



A++A_+(^++ip_ 



COS 



(20) 



where for brevity we introduced the notations A± = 



Here we pause for a 



I 7r/4 and (y5± 

moment to interpret Eq. (jl9p . From the discussion below 
Eq. (|T0)) and from Eq. ([T5|) it is clear that if we consider 
the two branches of Eq. ([9]) as classical Hamiltonians for 
two different quasiparticles then S^{0,y^) gives (half of 
the) enclosed flux by the quasiparticle trajectories with 
the wall between two subsequent collisions and 6*^(0) is 
the deflection angle of the momentum between the col- 
lisions. Therefore Eq. ([TOl) means that we neglect the 
second and higher powers of the difference between the 
enclosed flux and momentum deflection. 

It is clear that Eq. ([20)1 defines a pair of quantiza- 
tion conditions: A+ + A_ + ip^ + ip^ = (2m + 1)tt and 
A+ + A_-(p+-(f_ = {2l + l)nwithl,m = 0,l,2,.... By 



introducing the angle /?± = arcsin 



\r± 



). 



where X — kl . 



is the guiding center coordinate and R± — Q±/{eB) is 
the radius of the cyclotron motion for the two quasipar- 
ticle branch, we finally arrive at the following two quan- 



tization conditions: 
'1 



Rl 



and 
1% 



;sin2/3+ + /3+ 



i sin 2/3_ + ! 



(/3+ -f /3_) ^ Anm, 



TO = 0, 1, . . . rrimax 

(21a) 



1 



;sin2^+ + + - + 



2^ 



1 



■sin2/3_ + - + 



0. 



I — 0,1,... Imax 

(21b) 



These equations are the first important results of our pa- 
per. They do not lend themselves to a simple semiclassi- 
cal interpretation but a possible classical picture could be 
the following: the classical skipping orbits whose quanti- 
zation would be described by these equations consist of 
two segments, each of them having slightly different radii 
given by R± but the same guiding center coordinate. Be- 
sides the orbital motion, the quantization conditions also 
depend on the change of the phase of the spinor part of 
the wavefunction which is described by the (/3+ -I- /3_), 
(tt - 13+) + {tt- P^) terms in Eqs. ((2T|) . 

We note that analogous calculations to the ones out- 
lined above can be carried out if the dominant term in 
the SOI is the fc-linear Dresselhaus term. Therefore this 
approach can be relevant e.g. in the case of heterostruc- 
ture studied in Ref. iA^ 




FIG. 2: (color online) Comparison of the dispersion of the 
bands as a function of the wavenumber k obtained from tight- 
binding numerical calculations (dots), and using the theoreti- 
cal prediction given by Eqs. ^21\ (solid lines) for magnetic field 
B — 0.15T and ksols ~ 0.256 (other parameters are given in 
the text). The dashed line at the leftmost band shows the 
result of Eq. ((22]) for m = 0. 



To see the accuracy of the semiclassical quantization 
we have performed numerical calculations for the dis- 
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persion of the edge states using the tight-binding ver- 
sion of the Hamihonian ([IJ [see e.g. Ref. |^ for the ex- 
plicit form of the tight-binding Hamiltonian] . The re- 
suhs for ksols — 0.256 at a relatively weak magnetic 
field of ,B = 0.15T and using typical parameters of 
2DEG in InSb quantum well at higher electron densities^ ^ 
(to* = 0.021 TOe, where me is the bare electron mass and 
aR = 1.4 * 10~^^ eVm) are shown in Fig. [21 As one can 
see Eqs. ((2T|) describe quite well the dispersion of the 
subbands, even at low energies, except for the k values 
where the guiding center X « i.e. in the transi- 
tion region to the bulk Landau levels. We have found 
that although Eqs. (|2T|) have solutions even for to = 
and I = 0, the leftmost band in Fig. [21 which is related 
to the zeroth Landau level, is poorly approximated by 
any of the to = or Z = curves that can be obtained 
from Eq. (|21ap or Eq. (|21bp . respectively. Nevertheless, 
the approximation works quite well for rest of the sub- 
bands, i.e. for TO,Z > 1. For a stronger magnetic field 
of S = 0.6T {k.olB = 0.128) shown in Fig. [1 the ap- 
proximation for the TO = 1 and I = 1 bands deteriorated 
as well, while higher subbands are still well described by 
Eq. (HID. 




^ \ 

-4 -2 2 4 6 



FIG. 3: (color online) Comparison of the dispersion of the 
bands as a function of the wavenumber k obtained from tight- 
binding numerical calculations (dots), and using the theoret- 
ical prediction given by Eqs. \21} (solid lines) for magnetic 
field B = 0.6T and kaols ~ 0.128 (other parameters are given 
in the text). Dashed lines at the leftmost bands show the 
result of Eqs. ^ for m = 0, 1. 

It is interesting to note that we find that the bands 
for <C |X| < R± are quite well approximated by the 
following simple formulae: 

f-sin2/3±+/?± + -j -7re(X) = 27r(TO + 3/4) 

(22) 

where Q{x) = 1 for a; > and zero otherwise. Note that 
in contrast to Eqs. (|^T|) these equations give a semiclassi- 




FIG. 4; (color online) Explanation of the phase shift —tt for 
X > in Eq. (|22p . Solid (blue) lines show the bent quasipar- 
ticle trajectories. The direction of the momentum at the start 
of a full period and at the subsequent collision with the wall is 
indicated by arrows. Also shown are the angles that the mo- 
mentum encloses with the wall. The phase shift contribution 
coming from the spinor part of the semiclassical wavefunction 
is ^(9out ~ Oin + 27) in (a) and ^{9out - Oin - 27) in (b). 

cal quantization of the orbital motion of two independent 
systems whose classical motion is described by the Hamil- 
tonians given by the left-hand side of Eq. ([9]). The 
spinor nature of the quasiparticles enters the quantiza- 
tion only through a — tt phase shift for large positive X 
values. The origin of this phase shift can be understood 
by looking at Fig. [4l For —R± < X <C the phase con- 
tribution coming from the ± exp(i6'(y)/2) factors of the 
wavefunction (sec Eq. (jlip is zero over one full period of 
the classical motion [ Fig. [HJa)]. This happens because 
the phase change accumulated during the orbital motion 
Sout — Qi = —27 is canceled by the phase change 27 upon 
refiection when the sign of the 11° is negated. However, as 
it is explained in Fig. IHJb), for orbits with Q ^ X < R± 
the total phase change is —2n/2 — — tt. We find that the 
leftmost band in Fig. [21 which is related to the zeroth bulk 
Landau level (see Eq. [SJ can be quite well approximated 
by the — branch of Eq. [22lfor m = 0. For stronger mag- 
netic fields, such as shown in Fig. [31 the left-most bands 
are better approximated by the quantization (|22|) than 
by Eqs. (|^T|) . We note that for X < the quantization 
in Eq. (|22l) basically corresponds to the "longitudinal SO 
approximation" studied in Ref. [TtI 

Another interesting and important comparison of the 
quantization shown in Eqs. (12ip can be made to the 
closely related results of Ref. ^i, where the authors used 
a different semiclassical approaches to describe the edge 
states in the presence of RSO. Firstly, we find the same 
result for the fc = energy levels as in Ref. H: 

Era{k = 0) = hLOcin - {k^Je f) (23) 

and a comparison with the numerical calculations show 
that it is a good approximation of the exact results. How- 
ever, we obtained a pair of quantization conditions, not 
just one [see Eq. (27) in Ref. Q. Furthermore, for the 
parameter range we consider our results give a good ap- 
proximation of the numerically calculated bands not only 
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close to A; = but for the whole dispersion relation. Fi- 
nally, an important difference in the semiclassical inter- 
pretation of the skipping orbits is the following: in the 
classical picture put forward in Rcf. 8 the skipping orbits 
consist of two different type of arcs, having radii R± and 
guiding center coordinates X±. Moreover, the guiding 
center changes upon each reflection at the wall (see Fig. 6 
in Ref. §)■ In contrast, our approach tells that the guid- 
ing center remains the same throughout the motion. We 
think that this is physically plausible because the guiding 
center is a constant of motion. It is instead the reflec- 
tion angle that slightly changes at each collision with the 
wall as a consequence of having two Fermi surfaces with 
different radii. 



IV. MAGNETIC FOCUSING 



Having obtained the quantization condition for edge 
states in Eqs. ((2T|) . the calculation of the focusing mag- 
netic flelds B focus goes along the lines of the discus- 
sion of Ref. 0- We expect that the ballistic trans- 
port in a mesoscopic wire in the magnetic field regime 
where the cyclotron diameter is smaller then the wire 
width can be understood in terms of the edge states 
described in Section |llT] because they are the propagat- 
ing modes of this problem. If the injector is narrow, 
i.e. it is only a few Fermi wavelength wide, one can 
assume that it excites these modes coherently. There- 
fore, as long as the distance between the injector and 
the collector is smaller than the mean free path (and 
phase coherence lenght), the interference of the edge 
states can be important. Since the total wave func- 
tion of the system can be written as a sum of all popu- 
lated edge states |fcm), at a given Fermi energy, i.e. 

- EiCi$(2/,fci)e^'^'" + E„,C™l'(y,fc™)e^'='"- (here 
Ci, Cm are normalization constants), interference along 
the conflnement potential is determined by the phase fac- 
tors exp{ikmx), exp{ikix). Here the wave numbers km, h 
are determined by requiring that they satisfy Eqs. (|21ap 
and (j21bp respectively, for a given Fermi energy Ep and 
quantum numbers m and I. As in Ref. we assume 
that the current at the collector is determined by the un- 
perturbed probability density and therefore the focusing 
peaks are the results of the constructive interference of 
edges states with ki^m ~ at a distance x = L from the 
injector. This corresponds to the assumption that only 
electrons injected close to perpendicular direction can be 
focused onto the collector. It is convenient to introduce 
the following notation: R = i?+ — Ep = Ep/huc 
and we denote by kp — Q±{Ep)/h the radii of the two 
Fermi surfaces (circles) in the wavenumber space. Then 
for k kp we have 1 ^ fi± « k/k^ and expanding 
Eq. (|21ap in this small parameter we flnd that in lowest 



order 

h — 



4i? 



— V f^^™ - T;[kFR+ + kpR-] 

1 + 



Similarly, expanding Eq. (j21bp we obtain 



(24a) 



AR ^ 



27r(2/ + 1) - -[k+R+ + kpR 



(24b) 

In the semiclassical regime, where Ep ^ hujc and hence 
> 1 we can take (l + « (^1 - -J-^ « 1. There- 
fore the phase difference kiL—kmL at distance L from the 
injector between two edge states whose wavenumbers km 
and ki are given by Eq. p4a|) and Eq. (j24bp . respectively, 
reads 

kiL~ kmL ^2n(l- m)^ (25) 
2R 2R 



and we remind that R{B) = \/vio + '^m* Ep / eB. We see 
that if for a given magnetic field e.g. L — 2R{B) then the 
phase difference kiL — kmL will be an odd multiple of tt. 
This means that the spinor part of the wave function of 
these two edge states, ^{y, ki) and $(y, km) [see Eq. p^ ] 
will appear with opposite signs in the total wave function, 
which may lead to a near cancellation of \km) and |fc;). 

In general, whenever L is an odd multiple of 2_R, the 
phase difference will be an odd multiple of tt meaning 
that there may be a near cancellation between \ki) and 
\km)- On the other hand, the phase difference between 
edge states |fc„i) and \k.^i), both belonging to the same 
semiclassical quantization branch given in Eq. (I21ap is 



krr,L - k'L ^ 2ir 



L 



2R 



(26) 



and a similar expression can be derived for the phase dif- 
ference kiL — ki' L between edge states \ki) and {kf ) of the 
other quantization branch [Eq. (j21bl) ]. As long a.s E ^ 1, 
edge states given by the same quantization branch can 
interfere constructively at distances L = q2R{B), q = 
1,2,3, ... i.e. regardless of whether L is an even or odd 
multiple of the cyclotron diameter 2R{B). By construc- 
tive interference we mean that wave function of \km) and 
|fc^j) have the same global sign. (Note however that the 
-E ^ 1, i.e. Ep 3> fiJjJc condition gives an upper limit for 
the q values for which this reasoning is applicable). 

In contrast, if L is an even multiple of 2R (e.g. L = AR) 
we find from Eq. (|25|) that the phase difference kiL—k^L 
will be an integer multiple of 27r. This means that in this 
case not only edge states belonging to the same quanti- 
zation branch but also those belonging to different quan- 
tization branches can interfere constructively. Following 
the reasoning of Ref. 2? therefore we expect that there 
will be peaks in the focusing spectra for magnetic flelds 
B where L is even multiple of 2R(B) while we might not 
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see peaks at magnetic fields corresponding to L being odd 
multiple of 2R{B). This simple analysis then suggests 
that the foc using fields are given by integer multiples of 
B focus = 4:^^2771* Ef + P^so/eL. 

It is interesting to compare these predictions on the fo- 
cusing spectra to exact numerical transport calculations. 
Using the tight binding version^, of the Hamiltonian ((T|| , 
the transmission probability Tic{B) between the injector 
and collector was calculated by employing the Green's 
function technique of Ref. H^. The scattering region was 
of finite width W^^^ and was assumed to be perfectly bal- 
listic and infinitely long (see Fig. [T]). This means that 
the left and right ends of it act as drains which absorb 
any particles exiting to the left of right. The spin orbit 
coupling had a finite value in the scattering region but 
was set to zero in the injector and collector. To simulate 
the effect of finite temperatures we used a simple energy 
averaging procedure in the calculation of the transmis- 
sion curves: T^^{B) = J T„iB,E)(^-'2^§^'jdE where 

fo{E) was the Fermi function. The actual results shown 
in Fig. \5\ were calculated at T = IK temperature. As 
it can be expected, higher temperatures tend to smear 
the curves while at lower ones an additional fine struc- 
ture appears. Our calculations are very similar to those 
in Ref. except that we used slightly different contacts 
(see the inset of Fig. [T|). The contacts were always of a 
finite width (typically five-nine sites wide) and could in 
principle accommodate more than one (spin-degenerate) 
open channels. 

Figure. [SJa) shows the transmission Td a function 
of the magnetic field B. We used parameter values 
that approximately correspond to the measurements of 
Ref. |2l| on InSb quantum wells: electron density rig = 
3.25 * effective mass m* — 0.02 rrie and Rashba 

coefficient au = 1.4 * 10~^"'^eVm. It has been shown2& 
that the effective giromagnetic factor ge// of InSb is 
quite large and therefore the Zeeman spin splitting can 
be noticeable already at relatively weak magnetic fields. 
Therefore in our numerical calculation we took into ac- 
count the Zeeman term as well and assumed ge// = — 22. 
The distance L was 945nm and both contacts were tuned 
to accommodate one (spin-degenerate) open channel. We 
find that for these parameters the first focusing peak is 
split [see the right inset of Fig. ^a.)]. That for strong 
enough an the first peak is split was first noticed in 
Ref. |7j. The peak splitting in good approximation cor- 
responds to AB — « 10 mT. Comparing now the 



numerical result on the first focusing peak to our analyt- 
ical prediction, we see that the dip between the peaks is 
at the magnetic field value where the analytics predicts 
that destructive interference may take place for L ~ 2R 
(corresponding to B ^ 0.2T) but that the presence of the 
twin peaks at S_ = -0.196 T and B+ = -0.206 T is not 
explained by our approach. It appears that the most 
straightforward way to understand them is to assume 
spin-split cyclotron orbits, as in Refs.l7l lTTHT6l . Looking at 
the second focusing peak, one observes that it is located 
at L = 4i? (which happens for B = 0.4 T) quite accu- 
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FIG. 5: (color online) a) The transmission from the injec- 
tor to the collector as a function of the magnetic field. The 
right inset shows the close-up of the first focusing peak. The 
left inset shows that if one neglects the Zeeman interaction 
then the third peak is not split, b) the partial transmissions 
Tc+,i+ (sofid line) and Tc-,i- (dashed line), c) the partial 
transmissions Tc~^i+ (sofid line) and Tc+,i- (dashed line) - in 
this scale Tc-^i+ cannot be distinguished from Tc+^i-. Note 
the different vertical scale with respect to Fig. [S{b). 



rately, in accordance with our edge-state based theory. 
The fact that its amplitude is significantly larger than 
the amplitude of the first peaks seems to corroborate the 
theoretical prediction that in this case edge states be- 
longing to different quantization branches [Eq. (|21a[) and 
Eq. (j21bl) ] can constructively interfere with each other. 
Whether this enhancement of the amplitude could be ob- 
served in an actual experimental situation, however, de- 
pends on how specular the reflection at the confinement 
potential is. (Note that in the classical picture the second 
focusing peak correspond to trajectories which bounce 
off the boundary between the injector and collector once, 
see the dashed line in Fig. [T]) A small amount of diffuse 
scattering at the boundary may render the observation 
of this enhancement difficult. Finally, we find that close 
to L = QR [B = 0.597T), where our calculations predict 
that the wave functions of the edge states may cancel, the 
amplitude of the transmission is indeed small, but the fo- 
cusing peak at a slightly higher magnetic field is again 
not captured by our calculations. It seems that peaks 
which appear when L is odd multiple of 2R can not be 
described with the presented theoretical approach. The 
splitting of the third peak close to B = — 0.62T, reminis- 
cent of the splitting of the first one, is due to the Zeeman 
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interaction and not to the SOI. As we mentioned, because 
of the large ge/ / the Zeeman energy can be important at 
smaller magnetic fields than in e.g. GaAs. This is il- 
lustrated in the left inset of Fig. ^a) where we show the 
calculation for the third peak but without taking into ac- 
count the Zeeman term in the Hamiltonian. One can see 
that for the assumed strength of an = 1.4 * lO^^^eVm 
the third peak is not split in this case. 

In Fig. [Sfb) and (c) we show the partial transmissions 
assuming spin polarized injection/detection. Thus, e.g. 
Tc_|_^i_ refers to the transmission probability of electrons 
being injected in spin —1 eigenstate and collected in +1. 
In contrast to Ref. we chose as spin quantization axis 
the y direction (for the definition of the coordinates, see 
Fig. [U . The motivation to choose this axis comes from 
Ref. ,20, where it was shown that the average spin of the 
edge states (at least the low energy ones) pointed mainly 
in the direction perpendicular to the confinement poten- 
tial, i.e. along the y axis in our case. Comparing Fig.[SJb) 
and (c) one can observe that except for the first focus- 
ing peak, the spin-flip transmissions Tc+^i^, Tc^^i^ are 
always significantly smaller than Tc+^i+j Tc-a-- We also 
performed calculations (not shown here) where the in- 
jected electrons were polarized in the x direction, as in 
Ref.0, and we found that Tc+_i_, T'c-,i+ were, apart from 
the vicinity of the first peak, usually smaller in the case 
of ij polarized injection than for x polarized one. Further 



investigation of the average polarization of the spin of 
the edge states in the semiclassical limit and its effect on 
the partial transmissions in a focusing setup is left to a 
future work. 



V. SUMMARY 

In summary, we studied the role of edge states in trans- 
verse electron focusing setup for strong spin-orbit cou- 
pling. As a first step, employing a semiclassical approach 
we derived a good approximation for the dispersion re- 
lation of the edge states and briefly compared our re- 
sults to other approximation methods that can be found 
in the literature. We then studied the interference of 
the edge states as this is expected to have ramifications 
on the focusing spectrum. Comparison of our theoreti- 
cal results with numerical transport calculations suggests 
that certain properties of the focusing spectrum can be 
quite well understood in terms of the interference of the 
edge states. Nevertheless, the presented semiclassical ap- 
proach can not capture all the important characteristics 
of the transport calculations. Finally, we studied numer- 
ically the electron focusing when spin polarized injection 
was used. 
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